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ABSTRACT 

Frequency- domain methods are used to study the angles of arrival and 
departure for multivariable root loci. Explicit equations are obtained. 
For a special class of poles and zeros, some simpler equations that are 
generalizations of the single input-single output equations are presented. 
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1. Introduction 


The stud/ of the angles of arrival and departure for multivariable root 
loci began only recently. It has proceeded using two different approaches: 
the state-space approach used by Shaked, [3] which involves computationally 
arduous spectral decompositions of matrices; and the frequency-domain approach 
used by Postlethwaite [2], which involves the use of a Newton diagram to ob- 
tain a series approximation of the loci near the pole or zero of interest.. 

In this paper we take the frequency-domain methods of Postlethwaite and 
develop them further to obtain more detailed results. It is shown that, sub- 
ject to certain conditions, multivariable root loci depart from poles and 
arrive at zeros in Butterworth patterns whose orders come from the McMillan 
indices of the transfer-function matrix G(s) at the pole or zero in question. 
Explicit equations for these angles, requiring only the evaluation of poly- 
nomials at the pole or zero in question, are obtained, We also define a 
special class of higher-order poles and zeros for which much simpler 
equations may be used. These simpler equations turn out to be generali- 
zations of the single-input-single-output CSISO) root locus equations, 

The problem considered is the standard root locus set-up, in which a 
system represented by the transfer-function matrix G(s) is placed in a feed- 
back loop with a scalar gain k multiplying all channels. As k varies from 
zero to infinity, the closed-loop poles vary, and the plot of these vari- 
ations in the complex plane is the root locus. G(s) is an mxm rational 
matrix function of the complex variables, and is assumed to have full rank 
and be strictly proper. 
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Background 


The closed-loop poles are given by the solutions to 
AU.^det^I-GOo) . g " - « C .) /- 1 

- . . . + C-l) m det G(s) - 0 (1) 

wherd g * -1/k. The multivariable root loci are thus branches of the alge- 
braic function s(g)j however, for the purpose of determining the angles of 
arrival and departure we may regard the root loci as analytic functions of 
g. Multiplying (.D through by A^Cs), the least, common denominator of the non 
zero principal minors of all orders of G(s), we obtain 

*Cg,s) - A jn Cs)g m +A m _ 1 (sjg 111 ' 1 + ... + A 1 C.s)g+A o Cs)-0 (2) 

where the A^s) are all polynomials. The A.(s) will figure extensively 
in the results to follow. It may easily be shown ( see Yagle [5],p.,19) that, 
excepting single-point loci, the poles of G(s) are the zeros of A m Cs) and 
the finite zeros of G(.s) are the zeros of A (s) . 

The Newton Polygon 

The Newton polygon is a graphical device that can be used to find a 
series approximation to the function $(g,s) in the vicinity of a zero of the 
function. It is discussed in detail in Walker [4] j here we merely give in- 
structions for constructing the polygon associated to the function $(g,s) 
around s=0. 

Cl I Write A. (s) =b . s a i + (higher-order terms) for i*0.,.m 

a- x « 

(,2) Set up u and v axes and plot the m+1 points 
P i 3 C.u,v) = (i,a i ) , 1=0,... m. 

(.3) Join P Q to with a convex polygonal arc each of whose 
vertices is a P^ and such that no lies below the arc. 
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(4) For each pair of points P^, forming the endpoints of 

a segment, compute the slope z of the segment and solve the 

i Ik 

equation c J b.. + c b^*0 for c, If there is another point 

h i k 

P^ on this segment, solve instead the equation c b^+ c J b^+ c b^»0, 

(5) A series approximation to 4(g>s) in the vicinity of s*0 is 
then gafcs" 2 . 


After obtaining this series approximation to a branch of the root locus, 

* ■ 

we ray then compute its angle from Arg ^|- . This is the method used by 

• m 

Postlethwaite [2] . Although construction of a Newton polygon is not neces- 
sary to employ any of the results of this paper, the Newton polygon is the 
basis for all of the proofs given in the Appendix. 


3 . Main Results 

In this section we state the main result of this paper. 

Definition . Let the Smith-McMillan form of GC.s) be diag 
and let p by an nth-order pole. Let k^ be the largest integer 
such that C.s-p)kj [ d.. Cs) , j*l,...m. Then the (kJ are the 
structure indices associated with the pole p, Note that 

(l)zl k. ■ n >. (2) k 1 >k_>...>k . 

fal J l— 2— — m 

Structure indices for zeros are defined analogously. We may 
now state: 

Theorem 1 The root loci departing from an nth-oTder pole depart 
generically in Butterworth patterns whose orders are the non- zero 
structure indices of G(s) at the pole. For a pole p with non?zero 
structure indices [k^,k 7 ,,.,k ] the angles of departure are: 
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0 


depart, i,j * k. 


Arg 


•Cn^k.-k.-. . ,-k.) 

l_l| i_ A .(«) 

ds Cnt’kj-k^- . , ,-kjj 
d i.n-K r ,,,-K.- 1 j ^ 

ds ^ n ” k l~ ' ’ ’ m *i + l 


(s) 


+ j 360° 

k^ , j * 0,1, , . ,k^-l , i * l,2, 1lf r 
if and only if the following conditions are met? 

CD k x ^ k 2 ^ ... / k r 

, I'n'l w 

’ll 

w ji Cp) ••• w jj Cp) ] 


(2) det 


Cp) ... Wjj Cp) 


jt 0 , j ■ 1, , ,r 


|s-p 

(3) 

(4a) 

(4b) 


where W(s) A VCs)U(s) and U(.s) are unimodular matrices that trans- 
form G(s) in its Smith-McMillan form, i.e.,, 

/n. Cs) n (s}\ 

a t») ■ u(») di ^(3^r ••• rnr/ vf5) (5 ’ 

For the angles of arrival we have; 

Theorem 2 The root loci arriving at an nth-order zero arrive generi- 

cally in Butterworth patterns whose orders are the non- zero structure 

indices of G(s) at the zero. For a zero z with non-zero structure 

indices [k , ... k ., k ] the angles of arrival are: 
m-r m-i m 


® arrival, i,j “ k m _. 


j 560° 

k . 
m-x 


d^-VVi- 




. e C.n-k-k .-k .) i+1 

as m m-i m-i_ 


(s) 


d Cn "V 


j (n-k - . . , 
.ds m 


’ ” Vi+l 5 . 

C -T A. Cs) 

m-i+l J 1 


k .-1 , i*0,l, , , , 
m*i 


iS» 


+ 


# j*0,l, ,,, 


r C6) 


if and only if the following conditions are met: 


« k m * k „-l * ■ 

• * k m-r 

C.2), det 

1 

H 
• • H 

. « 1? M 
# 


V‘ J • 



«>>«• k mtI £ °- 
Proofs: See Appendix. 


(7a) 

j*m?l,m-2, ... m-r-1 (7b) 


It should be noted that whether or not the conditions (4a) and (7a) hold, 
the orders of the Butterworth patterns are given by the structure indices as 
long as conditions (4b) and (7b) or similar conditions (see Appendix) hold. 

It is not clear that these conditions hold generically; however, it has been 
shown (Byrnes and Stevens (1]) that a necessary and sufficient condition for 
the Butterworth pattern orders to be given by the structure indices is for 
certain matrices arising in the block-diagonalization of G(s) to have simple 
null structure. This imples that the conditions (4b) and (,7b) are generic. 


The following example illustrates the application of Theorems 1 and 2. 

Example : Compute the angles of departure for the root locus of 

s 2 + 8s + 17 s 3 + 10s 2 + 33s- + 34 

G(s) ^ s 3 + 9s 2 + 25s + 17 2s 4 + 21s 3 + 78s 2 + 117s 

s 4 + 4s‘ > + 8s + 8s- + 4 

+ 68 

It is straightforward to compute 

<Kg,s) - (s 6 + 6s 5 + 18s 4 + 32s 3 + 36s 2 + 24s + 8)g 2 

-(2s 6 + 25 s 5 + 125 s 4 + 325 s 3 + 493 s 2 + 420 s + 170)g 

+ (s 4 + 16s 3 + 98s 2 + 272 s + 289). = 0 
and to obtain (from A2(s))the third<-order poles -1 ± j . The structure 
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indices for these poles axe [2,1], and the condition (4b) is indeed 
satisfied. Thus loci depart from the pole -1+j in a second-order and a 
first-order Butterworth patterns with angles 


^depart , 1 ’ * Ar * 


+ n360® 


12s 5 + 125s 4 + 500s 3 + 975s 2 

+ 986s + 420 


120s 3 + 360s 2 + 432s 

m 

+ 192 

s»-l+j 


, n»0,l 


0 


61.8* , 241.8® 

s 4 + 16s 3 + 98s 2 + 98s 2 * 272s + 289 


depart , 2 


Arg 


12s 5 + 125s 4 + 500s 3 + 975s 2 + 986s + 420 


s« -1+j 


- 33.7°. 

By symmetry, the angles of departure from the pole -1-j are ?61,8*, 
118.2®, and -33.7°. 

4. Simple Poles and Zeros 

For a special class of poles and zeros the angles of arrival and depar- 
ture may be computed more readily using the following theorems instead of 
Theorems 1 and 2. 

Definition : An nth-order pole is said to be simple if its structure 

indices are [n,0, ... 0]. 

A similar definition applies for nth-order zeros. We now have: 

Theorem 5 Let the Laurent expansion of GCs) at an nth-OTder pole 
p be 

G(s) = - 1 t G _ + ... + G , + G„ + , , . (.8) 


(s-p) 


n -n 


(VP) -1 


Then, the pole p is simple if and only if tr G_ n is non-zero, and the 
angles of departure from p are 

8 dep« 'n* 1 * [’« i- 0 ’ 1 ’ 

- “ Arg (-(s-p) n tr G(s)| , j»0,l, ,.n-l, 


C.9) 
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Theorem 4 Let & be an nth-order zero of G(s), and let the 
Laurent expansion of G"*(s) at z be 

H -n + * 7ntT H -l * H o * (10) 

Then, if and only if tr H_ n is non-zero, the zero z is simple and 

the angles of arrival at z are 

® arrival * k ** [ tT M -n] * ^ '* "- 1 

• jArj [(s-r) n tr 0' 1 Cs)| s . I ]*^-j 2 - _ 
j»0, 1, . . . n-1. (11) 


Proofs ; See Appendix . 

For the case of m»2, we may simplify (11) by noting that in this 
case we have 


tr G’ 1 (s) 


_ tr GCs) 

" det G(s) 


(. 12 ) 


so that (11) becomes 

® arrival * n Arg ^(• s ’ z ^ 

j«0,l,,.n-l. Q2). 

Although Theorems 3 and 4 are easier to employ than Theorems 1 and 2, 
the most interesting thing about them is that they are striking generali- 
zations of the SISO root locus equations for computing the angles of arrival 
and departure. The only difference is that the trace of the multivariable 
transfer function matrix is substituted for the SISO scalar transfer function. 
Note that this generalization is observed only for simple poles and zeros, 
since in the SISO case all higher-order poles and zeros are simple, while in 
the multivariable case only some are. 

Further results are available for the case of first-order poles. It is 
easy to show (see [5],p.51) that the conditions for Theorems 1 and 3 are equi- 
valent, i.e.,tr G_ 1 is non-zero if and only if w^Cp) is non-zero. If this 


n tr G(‘s) 1 . j’360 B 

det G(s\ s ,J n 
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condition does nojt hold, we may use the following theorem: 

Theorem 5 Let the Laurent expansion of G(s) at a first-order 
pole p be (8) . Then if 

(1) tr G_ i ■ 0 , (2) tr (G.^) + 0 
the angles of departure from p is 

depart ' ClJ) 

2 

and the root locus branch departs as k (a border departure) , 

Proof : See [5], p.52. 

5. Conclusion The behavior of the anglos of arrival and departure for 
multivariable root loci has been studied from a frequency-domain point of 
view, and explicit equations for these angles have been obtained. Simpler 
equations are available for the case of "simple" higher-order poles and 
zeros, and it was noted that these equations are generalizations of the’ 

SISO root locus equations for angles of arrival and departure, All of these 
results depend on some genericity conditions, For fiTSt-order poles an 
equation for the angle of departure was given for a case wherein these 
conditions are violated. More work needs to be done in clarifying these 
assumptions, interpreting them, and obtaining equations for angles in cases 
where these assumptions are violated, 
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Appendix 

Proof of Theorem 1 We shall prove Theorem 1 by first: constructing the 


Newton polygon, and then using it to obtain the angle® of departure, A more 
detailed proof is given in [5] , 

From (1) and (5) we have 


A(g s) ■ det 


» det 


gI-G(s^-det|gI-UCs)diag(^i^) V(.s) 
gl -diag^^yj W(s) j 


CIS) 


so that the coefficients of (2) may be written as 

. — /principal minors of order 

h&*> * A. 0»)C-D 2-( /n t Cs)\ 

”' h " Vh of W(s] 


) 


A a c S )c-n 



corresponding^ 
principal J > h*l . . . m 
minor of W (s)/ 


(16) 


where we have used the Binet-Cauchy Theorem in (16) 

Now let p be an nth-order pole wi^h non- zero structure indices 
[k^kj, . . .k ] , and assume that p is not a single-point locus. Then we may 
write 

A m (.s) * (s-P) n A m (s) , A m (p) + 0 (17) 

Cs) ■ (s-p) k i 3 i (s), Cp) + 0, i»l...ra (18) 

and substituting these into (16) we get 

A . (s) - (s-i>) Cn ' k i"-'V A m Cs)C-l) h n lt s )"-n h Cs) 

"‘ h a 

( higher-order terms! 
k in s-p ) ' ' ' 1 

as) 
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det 


WuCs) ••• w a CW 


w hl C - s5 "* 


(recall - chat kj>k,> ... >k m 3 . Since n^O) and d^fo) are relatively prime, 
n^Cp) i 0 for i*i,* ,r. It should now be evident from (19) that the Newton 
polygon will take the form given in Fig. 1 if and only if the conditions (,4a, 

4b) are satisfied. Note, that if for example k.*k,, then (4b) is replaced 

* «* 

by the condition 


i^Cp) 

^ w n (!>) 


n,(P) 

+ 3jcp) "22^^° 


r-oi 


which we would still expect to be true in general. 


From the Newton polygon we may approximate the root loci depart- 
ing from p by the series 

k a c i (s-p) k i, i»l.,.r (21) 

where the coefficients c^ solve 

(,c.) 1-1 b iml +(= i ) i b i » 0 (22) 

and where b^^ is defined from 

a r.'x „ u (m -'i (n-k,-. , , -k, ) . (higher-order \ 

Vi (s) b i<>P) 1 i + \terms in s-pJ . (23) 

It should be evident from (21) that the departing loci are grouped into 
Butterworth patterns with orders (k^, i»l.,.r). The coefficients may be 
obtained from (23) by taking repeated derivatives, and this yields (3), 
proving Theorem 1, 

Proof of Theorem 2 This is essentially the same as the proof of Theorem 1. 
The only difference is that we now use i*m, m-1, ... m-r instead of 
i»l, 2, ... r. 
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Proof of Theorem 5 Taking the trace of and using (17) , we have 
, A^C.s} " " A ro ^ tr G ^* 

■ W V* J fe tr 0 V'") 


C24) 


and if tr G_ n is non-aero Theorem 3 follows immediately from Theorem 1, 

Proof of Theorem 4 The main diagonal elements of G’ l O) we the principal 
minors of order m-1 of G(s) divided by det G(s). So we have 


A faUk f s ')WP rinci P al minors of\ , n m-l 
a^sj %UUZ^ order m-1 of G(s) J u 


A m Cs)det G(s) £(main diagonal elements of G -1 ^)) C-l) m " 1 


» -A q Cs) tr G _1 C3) 


C2S) 


and the rest of the proof parallels the proof of Theorem 3. 
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